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Abstract 



We study the resummation of large logarithmic perturbative corrections to the partonic cross 
sections relevant for di-hadron production in hadronic collisions, H1H2 — > h\h%K y at high invariant 
mass of the produced hadron pair. These corrections arise near the threshold for the partonic 
reaction and are associated with soft-gluon emission. We perform the resummation to next- 
to-leading logarithmic accuracy, and show how to incorporate consistently cuts in rapidity and 
transverse momentum of the observed particles. We present numerical results for fixed-target and 
ISR regimes and find enhancements over the next-to-leading order cross section, which significantly 
improve the agreement between theoretical predictions and data. 



1 Introduction 



Cross sections for hadron production in hadronic collisions play an important role in QCD. They 
offer a variety of insights into strong interaction dynamics. At sufficiently large momentum transfer 
in the reaction, QCD perturbation theory can be used to derive predictions. The cross section 
may be factorized at leading power in the hard scale into convolutions of long-distance factors 
representing the structure of the initial hadrons and the fragmentation of the final-state partons 
into the observed hadrons, and parts that are short-distance and describe the hard interactions 
of the partons. If the parton distribution functions and fragmentation functions are known from 
other processes, especially deeply-inelastic scattering and e + e~ annihilation, hadron production 
in hadronic collisions directly tests the factorized perturbative-QCD approach and the relevance 
of higher orders in the perturbative expansion. 

Much emphasis in both theory and experiment has been on single-inclusive hadron production, 
HiH 2 — > hX [H [2], [3J, HI El d El - Here the large momentum transfer is provided by the high 
transverse momentum of the observed hadron. Of equal importance, albeit explored to a somewhat 
lesser extent, is di-hadron production, HiH 2 — > hih 2 X, when the pair is produced with large 
invariant mass M. In many ways, one may think of this process as a generalization of the Drell-Yan 
process to a completely hadronic situation, with the Drell-Yan lepton pair replaced by the hadron 
pair. The process is therefore particularly interesting for studying QCD dynamics, as we shall also 
see throughout this paper. Experimental data for di-hadron production as a function of pair mass 
are available from various fixed-target experiments P, QUI ITT] , as well as from the ISR [12J. On 
the theory side, next-to- leading order (NLO) calculations for this process are available [T3 | [TH [To] . 
They have been confronted with the available data sets, and it was found that overall agreement 
could only be achieved when rather small renormalization and factorizations scales were chosen. 
The NLO calculations in fact show very large scale dependence. If more natural scales are chosen, 
NLO theory significantly underpredicts the cross section data, as we shall also confirm below. 

In the present paper, we investigate the all-order resummation of large logarithmic corrections 
to the partonic cross sections. This is of considerable interest for the comparison between data 
and the NLO calculation just described. A related resummation for the single-inclusive hadron 
cross section [5] was found to lead to significant enhancements of the predicted cross section over 
NLO, in much better overall agreement with the available data in that case. 

At partonic threshold, when the initial partons have just enough energy to produce two partons 
with high invariant pair mass (which subsequently fragment into the observed hadron pair), the 
phase space available for gluon bremsstrahlung vanishes, resulting in large logarithmic corrections. 
To be more specific, if we consider the cross section as a function of the partonic pair mass m, 
the partonic threshold is reached when s = m 2 , that is, f = fn 2 /s = 1, where \/§ is the partonic 
center-of-mass system (c.m.s.) energy. The leading large contributions near threshold arise as 
a k s [ln 2A;_1 (l — f)/(l — f)] at the fcth order in perturbation theory, where a s is the strong coupling 
and the "plus" distribution will be defined below. Sufficiently close to threshold, the perturbative 
series will be useful only if such terms are taken into account to all orders in a s , which is what 
is achieved by threshold resummation [T6l \T7\ [T8| [T9] . Here we extend threshold resummation 
further, to cross sections involving cuts on individual hadron and the rapidity of the pair. 

We note that this behavior near threshold is very familiar from that in the Drell-Yan process, 
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if one thinks of rh as the invariant mass of the lepton pair. Hadron pair production is more 
complex in that gluon emission will occur not only from initial-state partons, but also from those 
in the final state. Furthermore, interference between soft emissions from the various external legs 
is sensitive to the color exchange in the hard scattering, which gives rise to a special additional 
contribution to the resummation formula, derived in [TTJ [20l El] . 

The larger f, the more dominant the threshold logarithms will be. Because of this and the rapid 
fall-off of the parton distributions and fragmentation functions with momentum fraction, threshold 
effects tend to become more and more relevant as the hadronic scaling variable r ee M 2 / S goes to 
one. This means that the fixed-target regime is the place where threshold resummation is expected 
to be particularly relevant and useful. We will indeed confirm this in our study. Nonetheless, 
because of the convolution form of the partonic cross sections and the parton distributions and 
fragmentation functions (see below), the threshold regime f — ► 1 plays an important role also at 
higher (collider) energies. Here one may, however, also have to incorporate higher-order terms 
that are subleading at partonic threshold. 

In Sec. |2] we provide the basic formulas for the di-hadron cross section as a function of pair 
mass at fixed order in perturbation theory, and display the role of the threshold region. Section [3] 
presents details of the threshold resummation for the cross section. In Sec. H]we give phenomeno- 
logical results, comparing the threshold resummed calculation to the available experimental data. 
Finally, we summarize our results in Sec. [51 The Appendices provide details of the NLO corrections 
to the perturbative cross section near threshold. 



2 Perturbative Cross Section and Partonic Threshold 

We are interested in the hadronic cross section for the production of two hadrons hyp, 

H 1 {P 1 ) + H 2 {P 2 ) - hy{Ky) + h 2 (K 2 ) + X , (1) 

with pair invariant mass 

M 2 = (Ky + K 2 ) 2 . (2) 

We will consider the cross section differential in the rapidities rjy, r\ 2 of the two produced hadrons, 
treated as massless, in the c.m.s. of the initial hadrons, or in their difference and average, 

Ar] = -(m-m), (3) 

V = ^(Vi + V2) ■ (4) 

We will later integrate over regions of rapidity corresponding to the relevant experimental coverage. 
For sufficiently large M 2 , the cross section for the process can be written in the factorized form 



M 4 dM2d/ ^ r)df j = zJ J dx a dx b dz c dz d f^(x a , fi Fi )f^ 2 (x b , ^ Fi ) z c D^(z c , fi F f)z d D^ 2 (z d , fi Ff ) 

m i da ab ^ cd / A n R fj. Fi fx Ff \ 

r, Ar),r],a s (/m r ,— ,— , , 5 
V m m m / 



drh 2 d/\rjdri 
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where r) is the average rapidity in the partonic c.m.s., which is related to r] by 

f) = fj- I In ( — ] . (6) 



2 \x b 

The quantity Arj is a difference of rapidities and hence boost invariant. It is important to note 
that the rapidities of the hadrons with light-like momenta K x and K 2 are the same as those 
of their light-like parent partons. The average and relative rapidities for the hadrons and their 
parent partons are also therefore the same, a feature that we will use below. Furthermore, in 
Eq. (jSJ) the f^' 2 are the parton distribution functions for partons a, b in hadrons H\ t 2 and D^f 
the fragmentation functions for partons c, d fragmenting into the observed hadrons h\^- The 
distribution functions are evaluated at the initial-state and final-state factorization scales [ipi an d 
HFfi respectively. denotes the renormalization scale. The da ab ^ cd / dfdfjdAi] are the partonic 
differential cross sections for the contributing partonic processes ab — > cdX', where X' denotes 
some additional unobserved partonic final state. The partonic momenta are given in terms of the 
hadronic ones by p a = x a P\, p b = x^Pi-, Pc = Ki/z c , pd = K 2 /zd- We introduce a set of variables, 
some of which have been used in Eq. (j3J): 

(7) 
(8) 
(9) 
(10) 

(11) 



s = 


(Pi + P2) 2 , 




T = 


M 2 
~S~ ' 




S = 


(x a Pi + x b P 2 f = 




rh 2 = 


\z c z d J 


M 2 

Z c Zd ' 


f = 


m 2 M 2 


T 


S XqX^ZqZ^S 


X a X\)ZcZd 



At the level of partonic scattering in the factorized cross section, Eq. ([51), the other relevant 
variables are the partonic c.m.s. energy and the invariant mass rh of the pair of partons that 
fragment into the observed di-hadron pair. We have written Eq. (jSJ) in such a way that the term 
in square brackets is a dimensionless function. Hence, it can be chosen to be a function of the 
dimensionless ratio m 2 / s = f and the ratio of rh to the factorization and renormalization scales, 
as well as the rapidities and the strong coupling. In the following, we will take all factorization 
scales to be equal to the renormalization scale for simplicity, that is, (jlr — Hfi — I^Ff = A*- We 
then write 



rh 4 da 



ab— >cd 



The variable f is of special interest for threshold resummation, because it is a measure of the 
phase space available for radiation at short distances. The limit f — > 1 corresponds to the partonic 
threshold, where the partonic hard scattering uses all available energy to produce the pair. This 
is kinematically similar to the Drell-Yan process, if one thinks of the hadron pair replaced by a 
lepton pair. The presence of fragmentation of course complicates the analysis somewhat, because 
only a fraction z c za of rh 2 is used for the invariant mass of the observed hadron pair. In the 
following it will in fact be convenient to also use the variable 

. rh 2 M 2 

S z c z d S 
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which is the ratio of the partonic m 2 to the overall c.m.s. invariant S and hence may be viewed 
as the "r-variable" at the level of produced partons when fragmentation has not yet been taken 
into account. This variable is close in spirit to the variable r = Q 2 /S in Drell-Yan. 

The partonic cross sections can be computed in QCD perturbation theory, where they are 
expanded as 

(14) 



'OC s 

^ab^cd = ( 

7T 



LO , a s NLO , 

u ab^cd "i ^ab^cd ~>~ 
71 



Here we have separated the overall power of 0(a 2 s ), which arises because the leading order (LO) 
partonic hard-scattering processes are the ordinary 2 — > 2 QCD scatterings. At LO, one has f = 1, 
and also the two partons are produced back-to-back in the partonic c.m.s., so that fj = 0. One 
can therefore write the LO term as 

u L ab t cd (f, A V , i))=6(l-f)6 (fj) wSL^Aij) , (15) 

where is a function of A?7 only. The second delta-function implies that fj = ^ln(x a /xb). 

At next-to-leading order (NLO), or overall 0(a^), one can have f ^ 1 and fj ^ 0. Near partonic 
threshold, f — > 1, however, the kinematics becomes "LO like". The average rapidity of the final- 
state partons, c and d (and therefore of the observed di-hadrons) is determined by the ratio x a /xb, 
up to corrections that vanish when the energy available for soft radiation is squeezed to zero. As 
noted in Ref. [22], in this limit the delta function that fixes the partonic pair rapidity fj becomes 
independent of soft radiation, and may be factored out of the phase space integral over the latter. 
This is true at all orders in perturbation theory. One has: 

uj^ cd (f,Ar},fj,a s (iJ,),iJ,/m) = 5 (fj) u a l b n t cd (f, Ar/, a 8 (fi), fi/rh) + uJ a f__ cd (f, Ar], fj, a a (ji),fi/m) , 

(16) 

where all singular behavior near threshold is contained in the functions ^XScd- Threshold re- 
summation addresses this singular part to all orders in the strong coupling. All remaining con- 
tributions, which are subleading near threshold, are collected in the "regular" functions u a e ^ cd . 
Specifically, for the NLO corrections, one finds the following structure: 

u$Srf(f,Ai^ )A i/ih) = 8(fj)[u a 1 £ cd (A V ,fi/m)8(l-T) 

+ u a ^°(f,A V ,f,,»/m), (17) 

where the singular part near threshold is represented by the functions cu^^ cd , u) a \'^} cd , uJ a \'^ cd , which 
are again functions of only Arj, up to scale dependence. The "plus" -distributions are defined by 

f(x)(g(x)) + dx= (f(x)-f(l))g(x)dx-f(l) g(x)dx . (18) 
xo J xq Jo 

Appendix A describes the derivation of the coefficients oj^'®^, oj^zI^ uj a \fl cd explicitly from a 
calculation of the NLO corrections near threshold. This will serve as a useful check on the 
correctness of the resummed formula, and also to determine certain matching coefficients. 

As suggested above, the structure given in Eq. (fT7|) is similar to that found for the Drell-Yan 
cross section at NLO. A difference is that in the inclusive Drell-Yan case one can integrate over all 
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At] to obtain a total cross section. This integration is finite because the LO process in Drell-Yan 
is the s-channel reaction qq — > £ + £~. In the case of di-hadrons, the LO QCD processes also have 
t as well as w-channel contributions, which cause the integral over Arj to diverge when the two 
hadrons are produced back-to-back with large mass, but each parallel or anti-parallel to the initial 
beams. As a result, one will always need to consider only a finite range in At]. This is, of course, 
not a problem as this is anyway also done in experiment. It does, however, require a slightly more 
elaborate analysis for threshold resummation, which we review below. 



3 Threshold Resummation for Di-hadron Pairs 



3.1 Hard Scales and Transforms 

The resummation of the logarithmic corrections is organized in Mellin- iV moment space |16j . 
In moment space, the partonic cross sections absorb logarithmic corrections associated with the 
emission of soft and collinear gluons to all orders. Employing appropriate moments, which we 
will identify shortly, we will see that the convolutions among the different nonperturbative and 
perturbative regions in the hadronic cross section decouple. 

In terms of the dimensionless hard-scattering function introduced in Eq. (|12l) the hadronic 
cross section in Eq. (J5J) becomes 

da H lH2 ->h lh2 X 



M 4 dM2dAvd - = E/ Q dx a dx b dz c dz d f^(x a )f^(x b )z c D h c ^z c )z d D h d %z d ) 

x uj ab ^ cd (f, Ar?,i7,a s (ju), , (19) 

where for simplicity we have dropped the scale dependence of the parton distributions and fragmen- 
tation functions. At lowest order, when the hard-scattering function u ab _ cd is given by Eq. f|T5]) . 
the cross section is found to factorize under "double" moments [23l [51], a Mellin moment with 
respect to r = M 2 / S and a Fourier moment in fj = fj + | hi{x a /x b )\ 

/oo nl do H\H 2 -^h\h 2 X 

■oo Jo M dM 2 dA V drj 

= E f^( N + 1 + W2)/f 2 (iV + 1 - iv/2)&*{N + 2)D h d \N + 2) 

abed 

x l~Jfj^ J\tt n -H{1-t) 5{fj) (^) 2 ^l c ,(Ar ? ) , (20) 

where the Mellin moments of the parton distributions or fragmentation functions are defined in 
the usual way, for example 

/f(JV)= [ 1 x N - 1 f a H (x)dx. (21) 



o 



We note that instead of a combined Mellin and Fourier transform one may equivalently use a 
suitable double- Mellin transform The last two integrals in Eq. (l2~U|) give the combined Mellin 
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and Fourier moment of the LO partonic cross section. Because of the two delta-functions, they 
are trivial and just yield the N and v independent result {a s /n) 2 uJ a Q b ^ cd {Ar]). One might expect 
that this generalizes to higher orders, so that the double moments 

I"" drje ivf > I df^Ua^^Ari^a^),^) (22) 

o — oo */ 

would appear times moments of fragmentation functions. However, this is impeded by the presence 
of the renormalization/factorization scale \i which must necessarily enter in a ratio with m = 
M / y/z c Zd. As a result of this dependence on z c and Zd, the moments D^(N + 2), D d 2 (N + 2) of 
the fragmentation functions will no longer be generated, and the factorized cross section does not 
separate into a product under moments. Physically, this is a reflection of the mismatch between the 
observed scale, the di-hadron mass M, and the unobserved threshold scale at the hard scattering, 
rh. Threshold logarithms appear when § approaches the latter scale, not the former. This implies 
that at fixed M there is actually a range of hard-scattering partonic thresholds, extending all 
the way from M at the lower end to v^S* at the upper. This situation is to be contrasted to 
the Drell-Yan process or to di-jet production at fixed masses, where the underlying hard scale is 
defined directly by the observable. 

We will deal with the presence of this range of hard scales rh by carrying out threshold resum- 
mation at fixed rh as well as at fixed factorization/renormalization scale. For this purpose, we 
rewrite the cross section (Tl9|) in a form that isolates the fragmentation functions: 



M * dM 2 dAr]dfj = ^ J dZc dZd Zc jD ^ c ' ^ Zd D d 2 ( z ^ /*) ^H^cd \t', Ar], fj, a s (/x), ^ 

(23) 

where again r' = rh 2 / S = rx a Xb and Qn 1 H 2 ^cd is given by the convolution of the parton distribu- 
tion functions and u a b^ c d- 

r 1 

n Hl H 2 ^cd (V, Ar],fj, a s (fi), — ) =y~* / dxadxbfa 1 {x a , fi) f b Hi (x b , fi) uJ a b^cd ( t, Arj, fj, a a (fi), — ) 

ab 

(24) 

with f] = f] — ^\n(x a /xb) as before. At fixed final-state partonic mass rh, the function fl^H^cd 
now has the desired factorization property under Fourier and Mellin transforms: 

[°° dfje ivf > I dr> {r') N ^n HlH2 ^ cd (r',Ar],f],a s ^),^) (25) 
J-oo Jo v mJ 

= f* 1 ^ + 1 + ZZ7 / 2 ' v)fb H2 ( N + 1 - ZZ/ A y) Uab-^cd (N, v, A V , a.{p), 



rh 



where 



uj a b->cd (N, v, Arj, a s (/i), 4-) = / df)e lur > dff N 1 u ab ^ cd (f, A77, fj, a s (fi), 4- ) • (26) 
\ mJ J_ 00 J v mJ 

Through Eqs. (T23T) — (T26|) we have formulated the hadronic cross section in a way that involves 
moment-space expressions for the partonic hard-scattering functions, which may be resummed. 
Because the final-state fractions Zi equal unity at partonic threshold, the scale rh in the short- 
distance function may be identified here with the final-state partonic invariant mass, up to cor- 
rections that are suppressed by powers of N. For the singular, resummed short-distance function 
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we therefore do not encounter the problem with the moments discussed above in connection with 
Eq. 



3.2 Resummation at Next-to-Leading Logarithm 

As we saw in Eq. ( Fl6l) . the singular parts of the partonic cross sections near threshold enter with 
5(f)). This gives for the corresponding moment-space expression 

^™ d (N,A V ,a s (ri,^) = / 1 #*"- 1 a^(*,A^a.( A i),£) ■ (27) 

which is a function of N only, but not of the Fourier variable v. Dependence on the Fourier variable 
v then resides entirely in the parton distributions. It is this function, w^™^, that threshold 
resummation addresses, which is the reason for the use of the label "resum" from now on. 

The nature of singularities at partonic threshold is determined by the available phase space 
for radiation as f — > 1. Denoting by the combined momentum of all radiation, whether from 
the incoming partons a and b or the outgoing partons c and d, one has 

(Pa+ Pb ) 2 (Pa+Pb) 2 ~V^' 1 } 

where k$ is the energy of the soft radiation in the c.m.s of the initial partons. 

At partonic threshold, the cross section factorizes into "jet" functions associated with the two 
incoming and outgoing partons, in addition to an overall soft matrix, traced against the color ma- 
trix describing the hard scattering [TTJ [20] . Corrections to this factorized structure are suppressed 
by powers of 1 — f. The total cross section is a convolution in energy between these functions, 
which is factorized into a product by moments in f ~ exp[— N(l — f)], again with corrections 
suppressed by powers of (1 — f), or equivalently, powers of N. This result was demonstrated for 
jet cross sections in |20j, and the extension to observed hadrons in the final state was discussed 
in [261 12Z]- The resummed expression for the partonic hard-scattering function for the process 
ab -> cd then reads [TZl HSI 1201 12H : 

ssna A v, aM, 4-) = a^ +i ( a .w, 4-) Af +i ( a .M, 4-) 

V m/ V m/ V m/ 

I J ab^cd V m / 

xAf+ 2 U(»), ^)Ar(« s M^) . (29) 
V m/ V mJ 

We will now discuss each of the functions and give their expansions to next-to-leading logarithmic 
(NLL) accuracy. 

The A^ (i = a, b, c, d) represent the effects of soft-gluon radiation collinear to an initial or 
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final parton. Working in the MS scheme, one has [TBI [T7J [THl 1201 !2"Tj : 



In A? (aM, ^ 



rhJ 



l z N-l 



1-Z 



+ 



q 2 



(l-z) 2 m 2 to 



-Ma s (q 2 ))\nN--B i (a s (q 2 )) 



Here the functions A4 and 5, are perturbative series in a s , 

2 



A (a. 



TT V 7T 



+ ... 



and likewise for Bi. To NLL, one needs the coefficients [28] : 



A 



(i) 



a, 



2 



C; 



-2nbn 



67 
18 



Ti- 



-N f 



(30) 



(31) 



(32) 



where Nf is the number of flavors, and 



C„ — Ci 



N 2 C -1 

2N r 



3. c g 



C A = N C = 3 



11C A - 2N 



f 



12tt 



(33) 



The factors A^ generate leading threshold enhancements, due to soft-collinear radiation. We note 
that our expression for the A^ differs by the iV-independent term proportional to B^' from that 
often used in studies of threshold resummation (see, for example, Refs. [THJ [29]). As was shown 
in [TTJ [20[ [21] , this term is part of the resummed expression and exponentiates. In fact, the second 
term on the right- hand- side of Eq. fl30l contains the large- N part of the moments of the diagonal 
quark and gluon splitting functions, matching the full leading power \i ^-dependence of the parton 
distributions and fragmentation functions in Eqs. f[23|) and f[25|) . We shall return to this point 
below. 



Each of the functions H ab ^ cd , <S/v,a&->cd> S a b^ c d in Eq. (120!) is a matrix in a space of color 
exchange operators [TTj, 20], and the trace is taken in this space. Note that this part is the only 
one in the resummed expression Eq. f[29|) that carries dependence on Arj. The H a i,^, cd are the 
hard-scattering functions. They are perturbative and have the expansion 

H ab ^ cd (A V , a .M, £) = fTjJL, (A,) + (A V , t.) + 0(a 2 ) . (34) 

The LO (i.e. 0(a 2 s )) parts H a ® b \ cd are known [TTJ [20J [21], but the first-order corrections have 
not been derived yet. We shall return to this point shortly. The S a b^ cd are soft functions. They 
depend on iV only through the argument of the running coupling, which is set to /jl/N [T7] . and 
have the expansion 

S M (A^, a., |) = S a °l cd + ^S a Z cd (At,, ^) + 0(a 2 s ) . (35) 



S 



The iV-dependence of the soft function enters the resummed cross section at the level of next- 
to- next-to-leading logarithms. The LO terms S a ° b ^ cd may also be found in [T71 [201 121] • They are 
independent of Arj. 

The resummation of wide-angle soft gluons is contained in the S a b-> c d, which are exponentials 
and given in terms of soft anomalous dimensions, T a b-^ c d'- 



SN,ab->cd (At?, a s (//), — ) = Pexp - / — T ab ^ cd (A77, a s (q 2 )) 

V m) 2 J ^2 q 2 



(36) 



where V denotes path ordering and where iV = Ne lE with 7^ is the Euler constant. The soft 
anomalous dimension matrices start at 0(a s ), 



T ab ^ cd (Ati, a s ) = ^ (At)) + 0(a 2 s ) . (37) 



7T 



Their first-order terms are presented in [T71 EQl EQ EOj . 

Note that the Born cross sections are recovered by computing Tr {H^ S(°'}ab->cd, which is 
proportional to the function uJ a ° b \ cd (Ar]) introduced in Eq. ( fl5l) . It is instructive to consider the 
expansion of the trace part in Eq. (1291) to first order in a s . One finds [31] : 

Tr {hS ] n SS n \ = Tr{H^S^} ab ^ cd + - Tr {- [H® (T^ + H^S^T^] In A 

I J ab—*cd 71 

+HWS® + H^S^} ab ^ cd + 0(a 2 s ) . (38) 
When combined with the first-order expansion of the factors A^ in Eq. (1291) . one obtains 

u>T%i(N,AT),a s (n),^) = Ti{H^S^} ab ^ cd (l + ^ £ A® [in 2 N + lnAln(/i 2 /m 2 

\ i=a,b,c,d 



Tr {- \H^(T^SW + H <®S®rM] In TV 

7T 

+H^S^ + H^S^} aM + 0(a 2 s ). (39) 

This expression can be compared to the results of the explicit NLO calculation near threshold 
given in Appendix A. This provides a cross-check on the terms that are logarithmic in N, that is, 
singular at threshold. From comparison to the part proportional to 8(1 — f) in the NLO expression, 
one will be able to read off the combination (H^S^ + S^) in Eq. (139|) . This is, of course, not 
sufficient to determine the full first-order matrices H^ 1 ' and which would be needed to fully 
evaluate the trace part in in Eq. (EI) to NLL. To derive and S {1 \ one would need to perform 
the NLO calculation near threshold in terms of a color decomposition [32], which is beyond the 
scope of this work. Instead, we use here an approximation that has been made in previous studies 
(see, for example, Ref. [5]), 

TrfHS^SsA ~ (l + ^in T r{^4s(%} (40) 

L ) ab^cd V 7T / L ) ab—*cd 

where 
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are referred to as "C-coefficients" . The coefficients we obtain for the various partonic channels 
are given in Appendix B. The approximation we have made becomes exact if only one color 
configuration contributes or if all eigenvalues of the soft anomalous dimension matrix are equal. 
By construction, it is also correct to first order in a s . 

We now turn to the explicit NLL expansions of the ingredients in the resummed partonic cross 
section. For the function A^ in Eq. ( 1301) one finds: 



InAf (a s (ji), 4) =hf\\) \nN + hf ] (x, a s {fi), 4-) + ln^ (a, a.(ji), 4-) , (42) 

where A = boa s (fi) IniV and the functions hf\ h\ 2 \\n(£i) are given by 

4 (1) 

h\ (X) = —\— (2A + ln(l-2A)) , 



2tt6 A 

2-nbn \ % 7ib n m z 



in 



+4\ 

4irb 3 



V m/ 7rbo \ 2 



bi 


4 2) 




TT&o 


2tt6 


ln(l 


}") 


ln(l 



m 2 



(43) 



We note that we have written Eq. (jl2"|) in a "non-standard" form that is actually somewhat more 
complex than necessary. For example, one can immediately see that the terms proportional to 
ln(l — 2 A) cancel between the functions hf^ and ln(£j), as they must because they were not 
present in the in Eq. fl30l) in the first place. The term proportional to ln(/x 2 /m 2 ) in ln(£j) 
is the expansion of the second term in Eq. ( 1301) . Its contribution involving does not carry 
logarithmic dependence on N and would normally be part of the "C-coefficients" discussed above. 
The term proportional to ln(l — 2A) in ln(£j) has been separated off the first term in Eq. (130]) . Our 
motivation to use this form of Eq. (f42|) is that the piece termed ln(£j) may be viewed as resulting 
from a large- N leading-order evolution of the corresponding parton distribution or fragmentation 
function between scales m/N and the factorization scale \ip (we remind the reader that we have 
set the factorization and renormalization scales equal and denoted them by As mentioned 
earlier, the factors (— 2A^ In N — B^) correspond to the moments of the flavor-diagonal splitting 
functions, P??, while the term in square brackets is a LO approximation to 

bo %a s (q>). (44) 

J Hp Q 



Therefore, it is natural to identify 

£t (A,a„(/J),4) l H (N^) ~ l H (N,m/N) , (45) 
V mJ 

that is, the exponential related to Si evolves the parton distributions from the factorization scale to 
the scale m/N, and likewise for the fragmentation functions. At the level of diagonal evolution, it 
makes of course no difference if ln(£j) is used to evolve the parton distributions or if it is just added 



10 



to the function h\ ' . However, as was discussed in [33], [33] , one can actually promote the diagonal 
evolution expressed by £ \ to the full singlet case by replacing the term (— 2Af' In N — B^) by the 
full matrix of the moments of the LO singlet splitting functions, P^ , so that £ itself becomes a 
matrix. Using this matrix in Eq. (H2l) instead of the diagonal £i, one takes into account terms that 
are suppressed as l/N or higher. In particular, one resums terms of the form ln 2fc_1 N/N to all 
orders in a s [33] • We will mostly stick to the ordinary resummation based on a diagonal evolution 
operator £i in this paper. However, as we shall show later in one example, the subleading terms 
taken into account by implementing the non-diagonal evolution in the parton distributions and 
fragmentation functions can actually be quite relevant in kinematic regimes where one is further 
away from threshold. Here we will only take the LO part of evolution into account, extension to 
NLO is possible and has been discussed in [55] . 

For a complete NLL resummation one also needs the expansion of the integral in Eq. (I36p . 
which leads to 

IhSa^c (At,, M/i), £) = H \~* X) rit c , (A17) • (46) 

As in [29], we perform the exponentiation of the matrix on the right-hand-side numerically, by 
iterating the exponential series to an adequately large order. 



3.3 Inverse of the Mellin and Fourier Transform and Matching 
Procedure 

As we have discussed in detail, the resummation is achieved in Mellin moment space. In order to 
obtain a resummed cross section in r space, one needs an inverse Mellin transform, accompanied by 
an inverse Fourier transform that reconstructs the dependence on fj. The Mellin inverse requires 
a prescription for dealing with the singularity in the perturbative strong coupling constant in 
Eqs. f|3Dl. fl36l) or in the NLL expansions, Eqs. fH2l) . fH3|) . We will use the Minimal Prescription 
developed in Ref. [35], which relies on use of the NLL expanded forms Eqs. fj42|) . fj43|) . and on 
choosing a Mellin contour in complex-iV space that lies to the left of the poles at A = 1/2 and 
A = 1 in the Mellin integrand. From Eqs. ( |25i) and ( 1261) . we find 

nSS-Krf Ar?, fj, a 8 {p), = due-** / (r')^ 

1 2 V mJ 2n Jc MP -ico 2ttz 

x £ fS 1 (N + 1 + iu/2, 2 (N +1 — iv /2, /*) SJES (N, is, Atj, a.(/*), £j , (47) 

ab 

where the Mellin contour is chosen so that 6 Q! s (/x^) In Cmp < 1/2, but all other poles in the 
integrand are as usual to the left of the contour. The result defined by the minimal prescription 
has the property that its perturbative expansion is an asymptotic series that has no factorial 
divergence and therefore no "built-in" power-like ambiguities [35]. Power corrections may then 
be added as phenomenologically required. For most of our discussion below, the resummed short- 
distance function Cj r ^ b s ^ d is specified directly by Eqs. ( H2l) and ( 143]) . When we refer to "full singlet 
evolution", however, we make the identification in Eq. (H5l) . and evolve the parton distributions 
and fragmentation functions to scale rh/N . In this case the exponential in uj r ^^ d is found from 
the h\ and h\ terms only in Eq. (|42|) . 
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We note that the parton distribution functions in moment space fall off with an inverse power 
of the Mellin moment, typically as 1/N 4 or faster. This helps very significantly to make the inverse 
Mellin integral in Eq. (j4Tj) numerically stable. In particular, the resulting functions ^H^-^cd are 
very well-behaved at high r'. This would be very different if one were to invert just the resummed 
partonic cross sections u^^d and attempt to convolute the result with the parton distributions. 
The good behavior of the ^nTfh^cd ma kes it straightforward numerically to insert them into 
Eq. (1231 . where they are convoluted with the fragmentation functions in terms of momentum 
fractions z at fixed rapidities. At this stage, it is straightforward to impose cuts in the transverse 
momenta and rapidities of the observed particles. This gives the final hadronic cross section 
M A d<y HlH2 ^ hlh ' 2X /dM 2 dAr]df]. We note that because of the presence of the Landau pole and the 
definition of the Mellin contour in the minimal prescription, the inverted ^Hj^-tcd nas support 
at t' > 1, where it is however decreasing exponentially with r'. The numerical contribution from 
this region is very small (less than 1%) for all of the kinematics relevant for phenomenology. 

When performing the resummation, one of course wants to make full use of the available 
fixed-order cross section, which in our case is NLO (C?(af)). Therefore, a matching to this cross 
section is appropriate, which may be achieved by expanding the resummed cross section to 0(a^), 
subtracting the expanded result from the resummed one, and adding the full NLO cross section. 
Schematically: 

^^.match | ^^.resum ^^.resum | | ^/^.NLO 

In this way, NLO is taken into account in full, and the soft-gluon contributions beyond NLO are 
resummed to NLL. Any double-counting of perturbative orders is avoided. 



+ da NLO . (4* 



4 Phenomenological Results 

We now compare our resummed calculations to experimental di-hadron production data given as 
functions of the pair mass, M. These are available from the fixed-target experiments NA24 [9] (pp 
scattering at beam energy E p = 300 GeV), E711 [10] (protons with E p = 800 GeV on Beryllium), 
and E706 pj] (pp and pBe with E p = 500 and 800 GeV), as well as from the ISR pp collider 
experiment CCOR p2] which produced data at y/S = 44.8 and 62.4 GeV. The data sets refer to a 
7r°7r°X final state, with the exception of E711, which measured the final states h + h + X, h~h~X, 
h + h~X with h summed over all possible hadron species. When presenting our results for this data 
set, we will follow [H] to consider for simplicity only the summed charged-hadron combination 
(h + + h~)(h + + h~)X. For this combination also the information on the fragmentation functions 
is more reliable than for individual charge states. 

In each of the experimental data sets, kinematic cuts have been applied. These are variously 
on the individual hadron transverse momenta px,i or rapidities r)i, or on variables that are defined 
from both hadrons, cos 8*, Y, p^ 11 . Here cos 8* is the mean of the cosines of the angles between the 
observed hadron directions and the closest beam directions, in a frame where the produced hadrons 
have equal and opposite longitudinal momenta, sinh rji = —pr,2 sinh r] 2 0, QUI dH [121 EH]- This 
system approximately coincides with the partonic c.m.s. In terms of the observed transverse 
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momenta and rapidity difference one has: 



cos 0* = -( m . —r , + Pr ' 2 - - A , J sinh(2A^) . (49) 

2 VPr,2+Pr,iCosh(2Ar7) p T>1 + p T>2 cosh^A??) J K u v 7 

Furthermore, y is the rapidity of the pion pair, 

where «; = Ifi + is the pair's four-momentum and where the second equality in terms of At/, fj 
and the hadron transverse momenta pr,i holds for LO kinematics as appropriate in the threshold 
regime. Finally, p^ air is the transverse momentum of the pion pair, 

pP air = |p T)1 + p T 2 | = \p T l - p T2 \ , (51) 

where again the second equality holds to LO. Thanks to our way of organizing the threshold 
resummed cross section, inclusion of cuts on any of these variables is straightforward. 

In all our calculations, we use the CTEQ6M5 set of parton distribution functions [3S], along 
with its associated value of the strong coupling constant. We furthermore for the most part 
use the "de Florian-Sassot-Stratmann" (DSS) fragmentation functions [37], but will also include 
comparisons to the results obtained for the most recent "Albino-Kniehl-Kramer" (AKK) set |38j . 
We note that one might argue that the use of NLO parton distribution functions and fragmenta- 
tion functions is not completely justified for obtaining resummed predictions, given that large- N 
resummation effects are typically not included in their extraction mostly from deeply-inelastic 
scattering (DIS) and e + e~ annihilation data, respectively. As was shown in Ref. [39] for the 
case of the Drell-Yan process, resummation effects in the parton distribution functions extracted 
from DIS appear to have a very modest impact, except when high momentum fractions and/or 
relatively low scales are probed, which is not the case for the data sets we are considering here. 
We expect the same to hold for the fragmentation functions. In fact, some large- N resummation 
effects have been included in the AKK analysis [3B], and comparisons to the results obtained for 
this set will therefore be interesting. 

We choose for our calculations the renormalization and factorization scales to be equal, and we 
give them the values M and 2M, in order to investigate the scale dependence of the results. One 
expects that a natural scale choice would be offered by the hard scale in the partonic scattering, 
which is 0{m). Because of the relation M = m^/z c z^ the scale M is actually significantly lower 
than m, typically by a factor 2. Our scale choices of M and 2M therefore roughly correspond to 
scales rh/2 and rh, and we refrain from using a scale lower than fi = M since this would correspond 
to a rather low scale at the partonic hard scattering. 

Figure [1] shows the comparison to the NA24 [9] data for pp — > 7T°n °I at = 23.7 GeV. 
The cuts employed by NA24 are |cos#*| < 0.4, average over \Y\ < 0.35, and p^ air < 1 GeV. 
We start by comparing the full NLO cross section to the first-order expansion of the resummed 
expression, that is, the last two terms in Eq. f|48l) . This will help to gauge to what extent the 
soft-gluon terms constitute the dominant part of the cross section, so that their resummation is 
reliable. It turns out that the two terms agree to a remarkable degree. The dashed lines in Fig. [1] 
show the NLO cross section for scales 2M (lower) and M (upper), while the crosses give the 
NLO expansion of the resummed cross section. Their difference actually never exceeds 1% for the 
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Figure 1: Comparison of the NLO (dashed) and resummed (solid (DSS) and dash-dotted (AKK)) 
calculations to the NA24 data JM/, for two different choices of the renormalization and factorization 
scales, /i = M (upper lines) and /i = 2M (lower lines). The crosses display the NLO 0(a s ) 
expansion of the resummed cross section. 
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Figure 2: Same as Fig. [7J but for charged-hadron production for pp scattering at y/S = 38.8 GeV 
and with cuts appropriate for comparison to El 11. The data are from J77y . 
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kinematics relevant for NA24. The solid and dash-dotted lines in the figure present the full, and 
matched, resummed results for the DSS and AKK fragmentation sets, respectively, including "C- 
coefficients" implemented as described in Sec. 13.21 (see Eq. (|4ip ). One can see that resummation 
leads to a very significant enhancement of the theoretical prediction. A very good description 
of the NA24 data [9] is obtained for both sets, much better than for the NLO calculation which 
falls short of the data unless rather low renormalization and factorization scales are used. Also 
the scale dependence of the calculated cross section is much reduced by resummation. We note 
that the resummed result for the AKK set shows a somewhat steeper M-dependence than that 
for the DSS set and lies lower at high M. This may in part be due to the fact that large-iV 
resummation effects were included in the AKK analysis of the e + e~ annihilation data, resulting 
probably in fragmentation functions that have an overall steeper ^-dependence. That said, given 
the still relatively large uncertainties of fragmentation functions overall, we also note that the 
different behavior of the AKK and DSS results might be just due to differing assumptions made 
in the respective analyses. 

We next turn to the cross section for charged- hadron production, pBe — > h^h^X, measured 
by E711 [10J at y/S = 38.8 GeV. We recall that we sum over the charges of the produced hadrons. 
The cuts applied by E711 were pr,i > 2 GeV, and average over —0.4 < \Y\ < 0.2. The cut on the 
individual hadron transverse momenta is, in fact, irrelevant for the values of M considered here. 
Furthermore, as stated in their Fig. 6 p33] for the pair mass distribution we apply p^ air < 2 GeV, and 
0.1 < | cos 6*\ < 0.25. Figure [2] shows the data and our results. As before, the agreement between 
NLO and the NLO expansion of the resummed calculation is excellent. Again, resummation leads 
to an increase of the predicted cross section and a reduction of scale dependence. Even though 
the resummed results agree with the data much better than the NLO ones for the scales we have 
chosen, they tend to lie somewhat above the data, in particular at the highest values of M. Keeping 
in mind the results for NA24, one may wonder if this might be in part related to the fragmentation 
functions for summed charged hadrons, which are probably slightly less well understood than those 
for pions, due to the contributions from the heavier kaons and, in particular, baryons. The trend 
for the resummed result to lie a bit high is, however, somewhat less pronounced for the AKK set 
which again produces results that are a bit steeper than the DSS ones. 

Figures [3] and 0] show the comparison of our results (for the DSS set) to the E706 data sets 
for neutral pion pair production in pp and pBe scattering at v^S* = 38.8 GeV (800 GeV beam 
energy), respectively. We do not take into account any nuclear effects for the Beryllium nucleus, 
except for the trivial isospin one. This has a very minor effect on the cross section, compared 
to pp. E706 used cuts fairly different from those applied in the data we have discussed so far. 
There were no explicit cuts on cos#*, p^ air or Y, but instead cuts pt,% > p'r 1 = 2.5 GeV and either 
-1.05 < rji < 0.55 (for the VS = 38.8 GeV data) or -0.8 < rji < 0.8 (for the v 7 ^ = 31.6 GeV 
data) on the transverse momenta and rapidities of the individual pions. The cut on transverse 
momentum, in particular, has a strong influence at the lower M: in a rough approximation, it 
leads to a kinematic limit M ~ 2px t i > 5 GeV, so that the cross section has to decrease very 
rapidly once one decreases M toward 5 GeV. This behavior is indeed seen in the figures. 

As in the previous cases, the NLO expansion of the resummed and the full NLO cross section 
agree extremely well, typically to better than 2%. For the two scales we have chosen, the NLO 
cross sections fall well short of the data. It was noted in [HI EE] that in order for NLO to match 
the data, very low scales of /i = 0.35M have to be chosen. The resummed cross section, on the 
other hand, has much reduced scale dependence and describes the data very well for the more 
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natural scales M and 2M, except at the lower M where the cut p™ 1 on the px t i becomes relevant. 
One observes that the data extend to lower M than the theoretical cross section, which basically 
cuts off at M — 5 GeV as discussed above. A new scale becomes relevant here, the difference 
\M — 2pj? t \. Higher order effects associated with this scale (which are different from the ones 
addressed by threshold resummation) and/or non-perturbative effects such as intrinsic transverse 
momenta [H] probably control the cross section here. It is also instructive to see that the cross 
section is very sensitive to the actual value of the cut on the py,i- In Fig. [5] we show the resummed 
results for scale \i = 2M for p T i > 2.5 GeV (as before) and pta > 2.2 GeV. One can see that with 
the lower cut the data are much better described. Experimental resolution effects might therefore 
have a significant influence on the comparison between data and theory here. 

In order to check consistency, E706 also presented their pBe data set at v^S" = 38.8 GeV 
when the E711 cuts were applied instead of the E706 default ones. These data are found in 
Figure [6] shows the comparison for this case. One can see the same trends as before. Clearly, the 
description of the data by the resummed calculation is excellent. For this set of cuts, the cross 
section is not forced to turn down by kinematics at the lower M, and theory and data agree well 
everywhere. Figures [7] and [8] show results corresponding to Figs. [3J HI but for the lower beam 
energy, 530 GeV, employed by E706 (y/S = 31.6 GeV). 

We finally turn to the data sets available at the highest energy, which are from the CCOR 
experiment at the ISR [12J. Two data set are available, at y/S = 44.8 GeV and 62.4 GeV. The 
cuts employed by CCOR were identical to those of NA24, | cos^*! < 0.4, average over \Y\ < 0.35, 
and p^ 1T < 1 GeV. Figure [9] shows our results at = 44.8 GeV. The resummed calculation again 
shows decreased scale dependence and describes the data much better than the NLO one. At the 
lower values of M, it does show a tendency to lie above the data. Barring any issue with the data 
(which appear to have a certain unexpected "shoulder" around M — 10 GeV or so), this might 
indicate that one gets too far from threshold for resummation to be very precise. On the other 
hand, the agreement between full NLO and the NLO expansion of the resummed cross section 
still remains very good, as can be seen from the figure. The trend for resummation to give results 
higher than the data becomes more pronounced at the higher energy, v^S* = 62.4 GeV, as Fig. [TUl 
shows, where we have used both the DSS and AKK sets of fragmentation functions. Although not 
easily seen from the figure, the NLO expansion of the resummed cross section starts to deviate 
more from the full NLO cross section than at the lower energies. At the lower M shown, it can 
be higher by up to 7%, which is still a relatively minor deviation, but could be indicative of the 
reason why the resummed result is high as well. 

Clearly, any deviation between the full NLO cross section and the NLO expansion of the 
resummed one is due to terms that are formally suppressed by an inverse power of the Mellin 
moment N near threshold. It is therefore interesting to explore the likely effects of such terms. 
This can be done by promoting the LO anomalous dimension in the evolution part in Eq. f|42l) 
from its diagonal form to the full one, as described in Sec. 13.21 

-24HnN-B\ l) ^PW> N , (52) 

which includes the sub leading terms 'ml/N and full singlet mixing. For simplicity, we perform this 
modification only for the lowest order part of evolution, as indicated in Eqs. fj4"3l and fl52l . The 
results obtained in this way are shown in Fig. [TTJ One can see that the resummed result obtained 
in this way indeed decreases significantly with respect to the one in Fig.[TU]which was based on the 
diagonal evolution only, and is much closer to the data. At the same time, the agreement between 
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Figure 3: Comparison of the NLO (dashed) and resummed (solid) calculations (for the DSS frag- 
mentation set) to the E706 pp data at yS = 38.8 GeV 171]/ , for two different choices of the 
renormalization and factorization scales, fi = M (upper lines) and \i = 2M (lower lines). The 
crosses display the NLO 0(a s ) expansion of the resummed cross section. 
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Figure 5: Resummed cross section for scale /i = 2M and p?,i > 2.2 GeV (dashed), compared to 
the one with pt,% > 2.5 GeV shown previously in Fig. \4\ (solid). 
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Figure 6: Comparison to E706 data with a different set of cuts, corresponding to the ones applied 
by E711. The data with these cuts are from [77] /. 
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the NLO cross section and the 0(a s ) expanded resummed result becomes as good as what we 
encountered in the fixed-target case. Figure [12] presents the corresponding result for the case of 
NA24. Comparison with Fig. [1] shows that the effect of the subleading terms is much smaller here, 
as expected from the fact that one is closer to threshold in the case of NA24. Nonetheless, the 
effects lead to a slight further improvement between the resummed calculation and the data. In 
particular, they give the theoretical result a somewhat flatter behavior, which follows the trend 
of the data more closely overall. While the implementation of subleading terms in this way will 
require further study, this appears to be a promising approach for extending the applicability of 
threshold resummation into regimes where one is relatively far away from threshold. 

That said, we remind the reader that already in the part that is leading near threshold we 
have made the approximation in Eq. (|41|) for our "C-coefficients" . This, too, will need to be 
improved in the future, by taking into account the full color structure of the hard scattering 
function beyond LO, as we discussed in Sec. 13.21 To give a somewhat extreme example of the 
effects generated by the C-coefficients, we have re-computed the resummed cross section for the 
case of CCOR at = 62.4 GeV, but leaving out all effects of the the coefficients beyond NLO. 
In other words, we leave out the C-coefficients in the first two terms on the right-hand-side of 
Eq. ( 1481) . keeping them of course in da NL °. This is likely not a good approximation of the beyond- 
NLO hard coefficients, because the C^_^ cd have ir 2 terms and logarithms in the renormalization 
scale fj, that are independent of the color channel and truly enter in the form given in Eq. (14T]) . 
Some of these are in fact even known to exponentiate [TTJ [201 EH HQ]. In any case, the 
result of this exercise is shown in Fig. [131 where it is also compared to our earlier calculation 
that included the C-coefficients in the way discussed in Sec. 13.21 One can see that there is a 
sizable numerical difference, and that the scale dependence of the resummed result without the 
beyond-NLO C-coefficients becomes significantly worse. 

We finally turn to the distribution in cos#*, defined in Eq. (1491) . for which most of the experi- 
ments mentioned above have presented data as well. In fact, the CCOR data [12] for this observable 
were instrumental in establishing the QCD hard-scattering nature of pp interactions [41] . From 
the point of view of threshold resummation, the distribution in cos 9* may appear somewhat less 
interesting than the pair mass one, since the threshold logarithms arise in 1 — f = 1 — m 2 /s, 
regardless of cos 6*. In addition, the cos^* distributions are presented as normalized distributions 
of the form 

da jd cos 9* 
da jd cos 0*|c OS 0*=o ' 

so that the main enhancement generated by threshold resummation is expected to cancel. Nonethe- 
less, as we have seen in Sec. I3.2[ the resummed expressions do contain additional dependence on 
An beyond that present in the Born cross sections, which will affect the cos 9* distribution at 
higher orders. This is visible from the soft part in Eq. fT46|) and also from the "C-coefficients" in 
Eq. (EH). Rather than going through an exhaustive comparison to all the available data, we just 
consider one example that is representative of the effects of threshold resummation on the cos 9* 
distribution. Figure [141 shows the normalized distribution for the E711 case, where we have again 
summed over all charge states of the produced hadrons. The dashed lines show the NLO result 
calculated again with the code of [13], for scales /i = 2M and fi = M. One can see that for these 
scales the NLO calculation is lower than the data for higher values of cos 9*. The dot-dashed lines 
in Fig. [T4] show the resummed results for scales /i = 2M and fi = M. These show a steeper rise 
with cos 9* and describe the data better than NLO for the scales shown. However, they still tend 
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Figure 9: Comparison of the NLO (dashed) and resummed (solid) calculations to the CCOR 
data [9] at \/~S = 44.8 GeV, for two different choices of the renormalization and factorization 
scales, /J, — M (upper lines) and fi = 2M (lower lines). The crosses display the NLO 0(a s ) 
expansion of the resummed cross section. 
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Figure 10: Same as Fig. 0, but for VS = 62.4 GeV. We also show the resummed result obtained 
for the AKK set of fragmentation functions. 
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Figure 11: As Fig. [70l but extending the diagonal evolution in the resummed formula to included 
subleading terms and singlet mixing, as shown in Eq. 15^) . We use the DSS set of fragmentation 
functions. 
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Figure 12: Same as Fig. [TTj, but for the case of NA24- 
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to lie below the data at higher values of cos 8*. As was suggested in [121 EH EES] , for the cos^* dis- 
tribution the hard scale in the partonic process will itself be a function of cos 9*, so that it is more 
natural to choose a factorization/renormalization scale that reflects this feature. We therefore 
present our resummed results also for scales /i = 2M* and [i = M*, where M* 2 = M 2 (l — cos 8*) 
which is proportional to the Mandelstam variable t in the partonic process. One observes that 
with these scale choices a very good description of the data is achieved. We note that in the 
NLO calculations presented in Refs. [HI [15] the scale was chosen proportional to the (average) 
transverse momenta of the produced hadrons, which for given M also depend on cos#*. This 
resulted in a satisfactory description of the data, when scales effectively a factor two smaller than 
our M* were used. Overall, the trend for the resummed cos 8* distribution to lie higher than NLO 
and be in better agreement with the data is found to be a generic feature that occurs as well for 
the cases of the other experiments. 



5 Conclusions 

We have investigated the effects of next-to-leading logarithmic threshold resummation on the cross 
section for di-hadron production in hadronic collisions, H1H2 — > hih^X, for a range of invariant 
masses of the produced hadron pair. We have developed techniques to implement the resumma- 
tion formalism at fixed rapidities for the produced hadrons and for all relevant experimental cuts. 
Extensions of these techniques to the level of next-to-next-lo-leading logarithms should be rela- 
tively straightforward in light of the close relation between the one- and two-loop soft anomalous 
dimension matrices [12]. 

For the fixed target and collider data studied here, the one-loop expansions of our resummed 
expressions approximate the corresponding exact one-loop cross sections excellently, to the level 
of a few percent and often less. In addition, with scales chosen to match the underlying hard 
scattering, the matched resummed cross sections typically explain the available data better than 
do NLO expressions at similar scales, with significantly reduced scale dependence. 

An important extension of these methods will be in the production and fragmentation of heavy 
quarks and in jet cross sections, where similar resummation methods are applicable. Given the 
reduction in scale dependence, this could provide an improved control over Standard Model tests 
and backgrounds in new physics searches. 
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Figure 13: Resummed results for the case of CCOR at = 62 A GeV. The solid lines show the 
results for scales M and 2M shown previously in Fig. Q/J, while the dashed ones were obtained by 
neglecting the contributions by the C^_^ cd coefficients beyond NLO. 




Figure 14: Normalized distribution in cos 6* (see for the case of charged-hadron production 

at E711. Dashed is NLO, while the dot-dashed and solid lines show resummed results. For the 
latter we have also used the scales /i = M* and /i = 2M* , where M* 2 = M 2 (l — cos#*). 
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Appendix A 



In this appendix we present some details for the calculation of the NLO partonic cross-sections 
near threshold. The virtual corrections have the 2 — » 2 kinematics of the Born terms and therefore 
fully contribute. They are proportional to 8(1 — r). The real-emission 2 — > 3 contributions require 
more effort. We consider the reaction a(pi) + 6(p2) — > c(k\) + d(k2) + e(ks), where partons d and 
e fragment into the observed pair of hadrons and have pair mass m 2 . It is convenient to work in 
the c.m.s. of the observed outgoing hadrons. We can then write the three-body phase space in 
4 — 2e dimensions as 

$ 3 = fA * n J ~) 25 T dTT-%1 - f) 1 - 2 * r dpp-%1 + p)- 2+2E 



(47r) 4 r(l -2e)\s J 

PTT 1*77 

x dip sin 1 "^ i[> d9 sin~ 2e 6 . (54) 



Here we define 



P = {pi- k 2 ) 2 /( P 2 - k 2 f = e" 2A " . (55) 



Near threshold, the integration variables are given in terms of the Mandelstam variables of the 
process as follows: 

(pi + P2) 2 = s , (k 2 + h) 2 = rh 2 = fs , 

(pi - h) 2 = - cos^) , (pa - ^) 2 = + cos^) , 

(Pi - ^ 2 ) 2 = ~j^- p = (P2 - h) 2 , (p, - fc 2 ) 2 = -j^— = (Pi - ^) 2 , 

( Ki + K2 ) = 1 + sinw cos 6 cos w 

v ' 2 V 1+P 1+P 

(fc 1 + fc 3 ) 2 = g(1 ~ r) f 1 - s in^cosg^- + cos^^^ ) . (56) 
2 V 1 + P 1 + P/ 



The phase space in Eq. (l5~4j) is used to integrate the squared 2^3 matrix elements \M. C 
For the latter one also assumes near-threshold kinematics. Since we want the partonic cross section 
at fixed f and A77, we only need to perform the last two integrations in Eq. f l54"l) . The basic integral 
for these is [43J 

' r ,„ sinl_2e V> sin" 2£ 6 

dip I dO 



(1 — cosipy(l — cos tfj cos x ~ sin-0 cos 0sinx) fc 

= 271 Y^T^ B ( l -z-^ l - e - k ) 2F1 (j, k,l~e, cos 2 1) , (57) 

where 2 -Pi is the Hypergeometric function. After integration over phase space and addition of the 
virtual corrections, infrared singularities cancel and only collinear singularities remain. These are 
removed by mass factorization, which we do in the MS scheme. Notice that since we are close 
to threshold only the diagonal splitting functions P-p contribute in this procedure. Combining 
all contributions, one arrives at the near-threshold structure of the partonic cross sections given 
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in Eq. (1171) . for each subprocess that is already present at LO. The final step is to take Mellin 
moments in f of the result, as described in Eq. fl26|) . This gives for the partonic cross sections to 
NLO: 



~thr,LO+NLO 



UJ 



ab^cd 



(N, A77, a s (p),p/m) 



CO 



(0) 



a s (n) 



71 



(1,0) 



->cd 



co^JAt}, p/m) - In N ^ cd ( Arp fi/m) + t (l n 2 N + C(2)) 



UJ 



(1.2) 



(At,) 



(5f 



where terms subleading in N have been neglected. The "C-coefficients" defined in Eq. (T4T1) are 
obtained from this as 



C$L* (Ar], p/m) 



(0) 

ab— >cd 



(Ar,) 



(59) 



Appendix B 



In this section we give the coefficients C a ]l^ cd for each subprocess contributing to the production of 
our di-hadron final state, resulting from the calculation outlined in Appendix A. In all expressions 
below, p is the renormalization scale. The dependence on the factorization scale is already included 
in the function Ei in Eq. ( j43l) . As before, we define p = e~ 2Av . 

qq' qq': 



We define: 
We then have: 



Q qq/ = 1 + (1+ P y. (60) 



C£L*(A,,,/*> = ^in£ + (^ + J§)tf, + (§-^)in»(i + ,) 



, 8 14 + 9p\ / 4 2 . , , 

+ ~o + hip+ -- + 777^ — ln(l + p)lnp 



3 QQqq' J \ 3 3Qqq 



P i \ 7vr 2 iV/ p 57V/ 

ln(l + p) + — — + — i- In 



3Qgg' 6Qgg' 3 1 + p 9 

8 T / p \ 3, , , 47tt 2 7 , . 

+-L 2 — — + - In 1 + p + + - . 61 



qq' qq': 



We have: 



^-L^(Ap,p/m) 



3 



7T 



;i + lnp) Inp + y 



Q 



qq' 



1 



2 Qqq' 

26 



ln(l + p) ln 



1 + P 



2Li 2 



1 + P 



(62) 



qq -> qq: 



We define: 



We then have: 



Q 



qq 



(l-p + p 2 )(3 + 5p + 3p 2 ) 
(l + p(l + p)) 



(63) 



2 ^ in ^ + 4 (l - p2)Li2 ( 



'7 



59p 5 



In p 



6 48Q OT 4Q g? 16 (3 + 5p + 3p 2 ) J 



(12p 2 + 3p-4) 



2g 



ln 2 (l + p) + 



qq 



'7 7 
M ^ - — (6 - 5p) 



lnp 

12Q qq 

53p-6 



37p - 71 + 



(17 - Bp) Q 



QQ 



3 4Q, 



99 



12 (3 + 5p + 3p 2 ) 



3 + 5p + 3p 2 
ln(l + p) lnp 



2 4Q gg 4 (3 + 5p + 3p 2 



/2-p 



ln(l +p) 



+ 



9 



2Q qq 3 (3 + 5p + 3p s 



,-) lnp-^Ml + p)-^ 



^ i + w - 



7T 



3g 



99 



7l7T 2 p 



16 V 144 (3 + 5p + 3p 2 ) ' 



(64) 



gg — > q'q': 



We define: 



We then have: 



C q ^ q , q , (Ar^p/m) 



Qg'g' = 1 + P 2 



27r6 ln4 + 7 ( ' 



3g, 



In 2 p 1 

1 12 1 



Qc 



ln 2 (l+p) 



7(1 + p) 



(65) 



6Q q 
5N f 



'Afop-ifi-*} ln(1 +p)lnp M 1+ s_m ln(1 + p) 



gg -> gg: 
We define: 



g« = 3 + p(l + p), 



l + 3p(l+p) 



(67) 
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We then have: 

CfiL* (A,, „/*) = 2,6, In £ + JV, (j + (1 + 2p) + j j In y t + f 



1 + P/V 2 0« 2 «« 3 / V 96Q« 2880!,? 36 

V-TO 61 + 180p 65 \ 2 / 13+124p 361 +972p 29. . 

, 7-p 35 + 71p ll\ , /61-64p 247 + 576p 22\ . 



, 8 + p 36 + 71p 7 \ , , , 5iV f 7 

+ — m + r/ + To ln 1 + p - -it + o ■ 68 

iegS 48Q® 12 / 9 2 



gg -> ##: 



We define: 
We then have 



G„ 5 = (l+p 2 )(4-p + V). (69) 



1 / 52p 2 p + 648 
+- 69+ ' 



) ln 2 (l + p) 



48 V G M - 4 - p + 4p 

15 . x2 9vr 2 (4-p) 191tt 2 14 

-p(l - p) 2 + — - 1 + — — - - . (70) 



4Ggq 16(4-p + 4p 2 ) 144 3 



We define: 



Q^2(l + p)+p 2 , 

g( 2 g )^9(l + p)+4p 2 . (71) 
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We then have: 



p 2 14 15(l + p)(2 + p)< 



2tt6 In - — + 



^Qqg Qqg 



2 /146 + 13p 3(109 + 13p) 241 



gg -» gg: 

We have: 



+ (1 + P) 



+ 



+ 



+ 



13 



15 



16Q<? 



13/0- 120 3(173 + 41p) 



24Q (1) 



to 



16Q (2) 



93 



86 + 89p 3(43 + 39p) 



6QS 



2Qqg 



17 
H 16 

27' 
16 



ln 2 p 
ln 2 (l + p) 



1 + P) 



89 

_ 3Qqg 



2)ln(l + p)+|-4 



31 27(p-3) 



&Qqg 



120 - 13p 3(155 + 23p) 31 . 
(TT~ W^ + 2A ' ln ' oln ( 1 + ' ) • 

YlQqg 8Q qg - ' 



c£Lqq (A77, n/m) = CjV gg (Ari,ii/m) 



99 -> 5 f 5 f : 
We define: 

= l+p(l + p) . 
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We then have: 



Cg> (An, li/m) = 2nb \n 



/i 



3P 2 (1+P)^ 
8 G lg 



m 

(3 + vr 2 ) - - 



5 3p 2 (l + p) 2 , 7i 2 p(l + p 2 )(l + pf 



9 



+ 



16 G lg 



3 /(l + p) J 



G 3 

S3 



+ 1 In 2 p 



^G 2 gg 



1 + P) 7^(! + P) +P -2(l + p) In l p 



G 



99 



Nt Nt 

(8 (1 + p 2 ) + hp) (1 + p) 2 ln(l + p) - — 4-(l + p) (5p 2 + 8(1 + p)) lnp 



24G 2 9 



24G 2 

99 



G 



-(l+p)-2-p ln 2 (l + p) 



.93 



8G lg 



(1 + p) 2p + 1 - — (1 + p) 2 lnpln(l + p) 



G 



99 



3tt 2 (1 + p) /3 3(1 + p) 3 , 3(2p + l)(l + p : 



4 G 2 



+ l 2 



^99 



+ 



AG 2 

^99 



lnpln(l + p) 



+ 



3(1 + p) 3 3(2p+l) 3(l + p) 3 



4 G 3 

93 



4G 2 

99 



2G 



99 



- l )\n 2 (l + P ) 



2G 



99 



tt 2 (1 + p) lbr 2 67 



2G 



.99 



6 
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